Fluxon Dynamics of a Long Josephson Junction 
with Two-gap Superconductors 
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We investigate the phase dynamics of a long Josephson junction (LJJ) with two-gap superconduc- 
tors. In this junction, two channels for tunneling between the adjacent superconductor (S) layers 
as well as one interband channel within each S layer are available for a Cooper pair. Due to the 
interplay between the conventional and interband Josephson effects, the LJJ can exhibit unusual 
phase dynamics. Accounting for excitation of a stable 27r-phase texture arising from the interband 
Josephson effect, we find that the critical current between the S layers may become both spatially 
and temporally modulated. The spatial critical current modulation behaves as either a potential 
well or barrier, depending on the symmetry of the superconducting order parameter, and modifies 
the Josephson vortex trajectories. We find that these changes in phase dynamics result in emission 
of electromagnetic waves as the Josephson vortex passes through the region of the 2-7r-phase tex- 
ture. We discuss the effects of this radiation emission on the current-voltage characteristics of the 
junction. 

PACS numbers: 74.20.De, 74.50.+r, 74.78. Fk, 85.25.Cp 



I. INTRODUCTION 

Recently discovered^ iron-based superconductors have 
renewed much interest in the Josephson tunnel junc- 
tions involving multi-gap superconductors*^ Both ex- 
perimental and theoretical studies indicate that sev- 
eral other superconductors, including MgB%^ NbSe 2 , 6 
heavy fermion UPt 3 £ organic^ (TMTSF) 2 X and K - 
BEDT, may have multiple pesudo-order parameters. As 
the appearance of a phase texture and unusual Abrikosov 
vortex properties^ in two-gap superconductors are con- 
sidered as manifestations of the multi-component order 
parameter, the tunneling property of a Josephson junc- 
tion with these superconductors may exhibit important 
differences from that with one-gap superconductors. This 
difference is attributed to the two tunneling channels for 
electrons in two-gap superconductor junctions due to the 
presence of two condensates. This internal freedom re- 
flects the number of electronic bands participating in su- 
perconductivity. Consequently, dynamics of the relative 
phase of the order parameters for the two-gap supercon- 
ductor tunnel junction differ from that of the one-gap 
superconductor junction. 

The presence of two tunneling channels indicates that 
there are two types of relative phase dynamics in a long 
Josephson junction (LJJ) with two-gap superconductors. 
These phase dynamics may be understood in terms of the 
interplay between the interband and conventional (intra- 
band) Josephson effects. The interband Josephson effect 
describes tunneling between two electronic bands in each 
superconductor (S) layer. On the other hand, the con- 
ventional Josephson effect describes tunneling between 
two adjacent S layers. These effects determine the dy- 
namics of the phase difference between the condensates 
within the same S layer and across two adjacent S layers, 
respectively. 

The relative phase of the two condensates is fixed in the 



ground state. This relative phase is locked to the value 
of and tt when the order parameter symmetry is S++ 

and 5-| , respectively. However, when the fluctuations 

about these phase-locked states are small, the interplay 
between the two Josephson effects can yield interesting 
phenomena. One such example is a collective excitation 3 
known as the Josephson-Leggett (JL) mode. This JL 
mode had been observed^ in MgB 2 by Bloomberg and 
coworkers, using Raman scattering. Also, theoretical 
studies of superconductor- insulator-superconductor junc- 
tions between one- and two-gap superconductors sug- 
gest that the ground state may violate^ the time rever- 
sal symmetry and that the phase dynamics of LJJ may 
depend on the symmetry of the superconducting order 



parameter 



12.13 



This hetero- Josephson junction may be 



fabricated— by using Nb (one-gap) and either MgB 2 or 
iron-based superconductors (two-gap). Ota and cowork- 
ers suggested that the gap symmetry can affect^ the 
Josephson current across the grain boundaries in poly- 
crystalline samples as well as the current-voltage (I-V) 
characteristics of the multi-gap intrinsic LJJ stacks. Re- 
cent theoretical studies of hetero- Josephson junction sug- 
gest that the phase dynamics of LJJ are affected by the 
JL mode. The effects due to availability of two tunnel- 
ing channels in LJJ may appear in measurable physi- 
cal quantities, including the drastic enhancement of the 
macroscopic quantum tunneling rate and the presence of 
an extra step structure,-^ 5 in addition to the conventional 
Shapiro steps, in the I-V characteristics. 

A deviation from the phase-locked state in LJJ may 
not be limited to a small amplitude. The experimental 
data from the magnetic response of a superconducting 
ring with two pseudo-order parameters indicate that a 
stable soliton-shaped phase difference between the two 
condensates (i.e., z-soliton) is attainable^ This observa- 
tion is consistent with a suggestion that the phase fluctu- 
ations can grow and produce a stable 27r-phase textured 
Excitation of i-soliton represents large phase fluctuations 



due to the interband Josephson effect and is taken as a 
hallmark of the multi-gap superconductors. 

Soliton states in two-gap superconductors had been 
explored by a number of authors. Kupulevaksky and 
coworkers examined^ this soliton state in mesoscopic 
thin-walled cylinders in external magnetic fields by us- 
ing the Ginzberg-Landau approach. Tanaka and cowork- 
ers suggested^ that interesting excitations, including a 
phase domain surrounded by the i-soliton wall, can arise 
in two dimensions (i.e., D=2) since an i-soliton may be 
considered as a D-l dimensional quantum phase disloca- 
tion. The i-soliton wall may carry a fractional flux quan- 
tum when one end of the soliton wall is terminated by 
the fractional vortex^ while the other end is attached to 
a sample edge. Also, a vortex-molecule may be formed 20 
when two fractional vortices, with a unit fluxoid quantum 
as the total magnetic flux, are connected by the i-soliton 
bond. These fractional vortices have been observed in a 
multi-layered superconductor— by using both magnetic 
force and scanning Hall probe microscopy. 

The i-solitons differ from the Josephson vortices (i.e., 
fluxons) since they do not carry magnetic flux and 
do not interact with either magnetic field or super- 
currents. However, an i-soliton may be formed and 
driven 22 by nonequilibrium charge density or by suffi- 
ciently strong superconducting currents. Gurevich and 
Vinokur suggested^ 2 , that spontanteous appearance of a 
soliton-like phase texture represents the breakdown of 
the phase-locked state. This breakdown can arise when 
the applied current density along the superconductor lay- 
ers exceeds the critical value. The phase fluctuations 
may appear as either an additional resonances in the AC 
Josephson effect or a static 27r-kink in the phase differ- 
ence. If the 27r-phase texture exists in each S layer, then 
this «-soliton may change the phase dynamics of the L JJ 
by inducing a critical current density modulation. 

Earlier studies on the effects of critical current mod- 
ulation indicate that both spatial and temporal depen- 
dence of Josephson current amplitude may be obtained 
by using experimental techniques such as Ohmic heat- 
ing, quasiparticle injection, and illumination with an in- 
tensity modulated beam of lights Also, a spatially peri- 
odic modulation of the critical current may be obtained 
from a periodic array of microresistors in the insulator 
(I) layer. 24 As the microresistors behave as pinning cen- 
ters for moving Josephson vortices (i.e., fluxons), the 
speed of the fluxon becomes modulated near each mi- 
croresistor. The effects of a small periodic critical current 
modulatio n 23 ' 25 " — created by an array of microresistors 
can yield a number of interesting properties, including 
emission of electromagnetic (EM) radiation and Joseph- 
son steps in the I-V characteristics! 23 ' 25 This suggests 
that the spatial and temporal periodic modulation of the 
critical current due to excitation of the i-solitons may 
exibit the same interesting property. 

Emission of EM radiation from an inhomogenous LJJ 
had been studied by using a number of different theoret- 
ical approaches, including the Green function perturba- 



tion technique^ 5 - inverse scattering perturbation theory 2 ^, 
and numerical simulation^ These studies reveal that a 
moving fluxon can radiate EM waves when its speed is 
larger than the critical value Vth- Also, the interference 
between the emitted EM waves can give rise to well- 
discernible steps in the I-V characteristics. This suggests 
that, if the critical current modulation is generated by 
large amplitude fluctuations of the phase difference, this 
may affect the junction property through the changes 
in the phase dynamics. These changes may be measured 
from the I-V curves. However, the effects of the interband 
Josephson current on the phase dynamics of a LJJ with 
two-gap superconductors have not yet been explored. 

In this paper, we consider a quasi-one dimensional LJJ 
(i.e. D=l) with the dimensions, compared to the Joseph- 
son length Xj, of L x ^> Xj and L y -C Xj (see Fig. 1). We 
assume that there are no fractional vortices and investi- 
gate the effects of large phase functuations (i.e., i-soliton) 
on the phase dynamics of the LJJ. We note that, here, 
a weak external magnetic field, applied parallel to the 
insulator layer, penetrates the junction in the form of 
fluxons. Before proceeding further, we outline the main 
result. We find that i) large fluctuations in the relative 
phase of the two condensates in each S layer via the in- 
terband Josephson effect may be described by the sine- 
Gordon equation, ii) The soliton-like excitation (i.e., 2n- 
phase texture) can generate both a spatial and temporal 
modulation of the critical current between the adjacent 
S layers, iii) The critical current modulation yields emis- 
sion of EM waves which radiate along the junction layer 
in a form of quasi-linear wave. The strength of the crit- 
ical modulation is characterized in the I-V curve as a 
discontinous structure. 

The outline of the remainder of the paper is as follows. 
In Sec. II, we describe the phase dynamics of a LJJ 
with two-gap superconductors by using a set of two sine- 
Gordon equations. In Sec. Ill, we derive the equation 
of motion for the relative phase of the two condensates 
in each S layer. In Sec. IV, we discuss the effects of 
the interband Josephson current on the phase dynamics 
by computing the radiation correction to the bare soli- 
ton solution. In Sec. V, we compute the trajectories 
of fluxon in the velocity-position phase plane and calcu- 
late the current- voltage characteristic curves. Finally, we 
summarize the result and conclude in Sec. VI. 



II. THEORETICAL MODEL FOR PHASE 
DYNAMICS 

In this section, we derive the equation of motion for 
the phase differences. Here, we describe the phase dy- 
namics of the LJJ by neglecting the dissipation effect, 
for simplicity, but this effect is included later in the cal- 
culation of the fluxon trajectories in Sec. V. Also, we 
neglect the boundary effect by considering a region away 
from the junction boundary in the x-direction (see Fig. 
1), where the Josephson critical current density J° in the 



absence of i-soliton is spatially not uniform^ We discuss 
the boundary effect of vanishing Josephson critical cur- 
rent on the I-V curves in Sec. V. 

We start with the following model Hamiltonian: H = 
Y)p Htb l + Ht- Here the Hamiltonian Htb i accounts 
for the two-gap superconductivity in the ^-th S layer, 
while the Hamiltonian "Ht accounts for electron tunnel- 
ing between the two adjacent S layers. First, we consider 
these Hamiltonian contributions in the absence of elec- 
tromagnetic fields. We write the two-gap Hamiltonian 
Htb, i as 




W-TB,, 



dr 



E 






nit 
Tin 



where e l describes the energy of electrons in the i-band 
(i = s, d) about the Fermi energy. For definiteness, we 
denote the two electronic bands involved in superconduc- 
tivity as s and d bands. The pairing interaction between 
electrons in the £-th S layer is described by the Hamilto- 
nian 'H p e air written as 



(1) FIG. 1. A LJJ with two-gap superconductor represented by 
two pseudo-order parameters *&Jexp(i# s ) and $o exp(i8 d ) is 
shown schematically. Here, L x and L y denote the dimensions 
in the x— and y— direction, respectively. J B is the bias cur- 
rent density and B is the externally applied magnetic field. 
d s and d denote the thickness of the superconductor (S) and 
insulator (I) layer, respectively. 
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V sd (c^cll e c d u c d u + h.c.) : (2) 



where c a e (c l a t ) denotes the operator which creates (de- 
stroys) an electron with spin a in the i-band, and V%j 
denotes the strength of pairing interaction between elec- 
trons in the i and j bands. On the other hand, the Hamil- 
tonian T-Lt describes the electron tunneling between the 
two adjacent S layers that are separated by the I layer. 
This Hamiltonian is given by 



Kt = J2 i T A<2 + h.c.) 



(3) 



»#:?> 



where T^- denotes the tunneling matrix element for an 
electron from the j to i band. 

Following Leggett^ we consider the subspace spanned 
by the BCS type function and introduce pairing operators 



&i — J\ Jt 



(4) 



to account for the two (i.e., s and d) condensates. Using 
the eigenvalues of the pairing operator ^>\ , we rewrite the 
Hamiltonian T-Ltb.i as 



h tb/ = nmf) + ti(\*i\ a ) - k s i*ii 2 - v M \* 



r/|2 



(5) 



where /Kl^l 2 ) corresponds to the kinetic energy of the 
electrons in the z-band. In the absence of contribution 
from the interband pairing (i.e., V sa - — 0), the electrons 
in the two bands (i.e., s and d) do not interact. The 
independent one-band gap equation for the parameter 
A 1 = Vu^ 1 may be obtained by minimizing Eq. ([5]). We 



may account for the phase effects of the two condensates 
by writing the complex pesudo-order parameter "J" as 



m 



K,ie 



i6\ 



(6) 



In this representation, the interband pairing interaction 
term in Eq. ([5]) depends explicitly on the relative phase 
of the two condensates (i.e., s and d) 



Xi 



(7) 



The interband pairing contribution describes the Joseph- 
son effect within each S layer due to tunneling of the 
condensates between the two bands. Similarly, the con- 
ditions for the energy extremum yield the two coupled 
gap equations of the form 



E^M 



(8) 



The coupled gap equations of Eq. ([8|) have two nontrivial 
solutions: \i = (i.e., Of = Of) and xe = ^ (i- e -> #f = 

9$ + ir) corresponding to the S++ and S_| symmetry, 

respectively. When V s d > 0, \e — is the stable solution. 
On the other hand, when V sc i < 0, xi — n i s t ne stable 
solution. 

In addition to the phase-locked state, soft modes as- 
sociated with fluctuations of xt-i representing a phase 
texture, may appear as a 27r-kink in xi- These soft 
modes can modify the phase dynamics of LJJ. For ex- 
ample, these modes can manifest as resonances in the 
AC Josephson effect when the two electronic bands are 
out of equilibrium. Also, they may appear as 27r-kink 
representing an i-soliton. 

We now investigate the effects of these soft modes by 
starting with the partition function Z, = Tr exp[— fil-L] for 
the LJJ with two-gap superconductors 29 in the presence 



of electromagnetic fields 

Z= J V[A°,A]V[^ l \D[6 l ]e- s[ ^°' el > A °- A] , (9) 

where /3 — 1/T and T denotes temperature. Here, we 
set h = c = kg = 1 for convenience. The scalar and 
vector potential in the ^-th layer are denoted by A° t and 
A( = (A x ,A V ,A z )t, respectively. The effective action 
S = Sg a p+Sfi e id~Tr G _1 is obtained by carrying out the 
imaginary-time functional integral over the Grassmann 
fields cj e and d a ^. The contributions from the psudo- 
order parameter and electromagnetic fields are denoted 
by S g ap and S field, respectively. The expression for S gap 
may be found in appendix and Sfi e M is included in the 
phase contribution, S p hase, to the action below. Also, we 
do not write the explicit expression of the inverse Green 
function G~ l here but note that it is a 8x8 matrix which 
consists of a 4x4 matrix for each two-gap superconductor 
layer. We extract the superconducting phase degree of 
freedom by performing the unitary transformation for the 
Green function and follow the standard procedure 3 -^ of 
retaining only the second-order tunneling processes aris- 
ing from the Josephson effect. The usual imaginary-time 
functional integral approach 3 ^ in this approximation al- 
lows us to obtain the phase contribution S p hase to the 
effective action S, which includes the phase terms, in an 
external magnetic field B applied in the y-direction (see 
Fig. 1) as 



Sphase 



P 



dr / dx Cr> 



(10) 



The effective Lagrangian density C p for the supercon- 
ducting phase s 31 ' 32 is given by 



C: 



8ire 2 



E 



i (de\ AO \ 2 i (de\ , X N 2 



y~] -r- cos Lp 13 + y^ mt _ er cos xe + Cem ■ (H) 



where e — 2e, and d s is the superconductor layer thick- 
ness. The phase difference of the superconducting order 
parameter in the magnetic field is yj IJ = 9\ — Q 2 — eA\ 2 , 

where A\ 2 = J. =1 A z (z)dz. The charge screening length 
Hi = yJXrpp/Airei is a constant related to the Thomas- 
Fermi screening length Xpr = \/ira /AkF, and the dielec- 
tric constant ti of the S layer. Here, a is the Bohr radius 
and Uf is the Fermi vector. The magnetic penetration 
depth of the S layer are denoted by \ — v/ m° j iirnie 2 
where m° is the mass of the i-band electron, 



' ^J^-^lTlVl ^ 



^ -^ j tanh -f + EU(El)[l - /(££)] 

k k 



El 



^ 



(12) 



A^, El = f3El, and f(E) is the Fermi func- 



tion. Here, we make a local approximation and consider 



/i and A as constants. 30 ' 31 The Josephson critical current 
J lJ betwen the electronic bands i and j in two adjacent 
S layers may be writte n 30 ' 31 as 



J' J 



2e 

'd^ 



*-2>-** 



k.k' 



%3 ElEl 




f(K>)} (13) 



[i-f(K)-f(El)] 



where u„ = 2nir/f3 with n = 0, ±1,±2, • • • is the Mat- 
subara frequency. The interband Josephson critical cur- 
rent J in ter = 2eV sd ^ s $> d / {V ss V d d - V 2 d ) between the two 
bands within the same S layer does not represent charge 
transfer, as in J y . The Lagrangian density for the elec- 
tromagnetic field Cem in the insulator layer is given by 



CEM = ^{e{El 2 f-{B\ a f] 



8- 



(14) 



where d is the I layer thickness, and e is the dielectric 
constant. The electric and magnetic field between two 
adjacent S layers (i.e., £ — 1 and 1 — 2) are defined as 



dA z 1 
El 2 = -^-- {Al -Al) 



1 

d { 



Bl 2 = -M X i- A 2) 



dx 



(15) 
(16) 



respectively. 

We now derive the equations of motion for the phase 
difference ip lJ by following the usual approach of mini- 
mizing the action S p hase of Eq. (flO)) . Here, noting that 
r = —it, we examine the phase dynamics of both ip ss and 
ip dd since we may write that ip sd = Of — 6 d + edAf 2 = 
ip ss + X2 and ip ds = 6 d - #1 + edA\ 2 = ip ss - X i- Us- 
ing the Euler-Lagrange equation for the scalar and vector 
potentials (i.e., A^ and A z ), we obtain the following two 
coupled equations of motion for the phase differences ip ss 
and ip dd : 



d 2 ip S! 



: 9_V 
dx 2 



2,„dd 



a> 8 2 if 

Cm ^9. + Ca 



dt 2 



dx 2 



2J S 



■ sin ip 



^ [sin(^ + X 2) + sin(^ ss - X i)](17) 
e 



Jinte 



■(sinxi -sinx2) = , 



and 



'/, 



d #V 



2,„dd 



+ - 



2J S 



" dt 2 

■sm<p dd + 



4 



d 2 if dd 

dx 2 



i2 . _ss 



c 



d 2 ip 



2,„ss 



dt 2 



+ C 



d (p 



dx 2 



J s 



sin(/ d - X 2)+sin(/ d + Xl )}l8) 



inter 



(sinxi -sinx2) = 



The coeffiecients in the above two equations of motion 
(for (p ss and ip dd ) are 



% 



4ire 2 £ 



2<r2 



and 



where 



Ce 






\47re£ s £ d 



d s d 

47t£ 2 K ? 

" d 



Kc 



4-7T 



1 



1 



+ k| 



(19) 



(20) 



(21) 



Here, £j denotes either the screening length (/i^) or the 
magnetic penetration depth (Aj) of the i-band, k° = e/27r 
and k£ = 1/27T. 

By noting that the phase difference y sd may be written 
as either yj srf = (p ss + \2 or (p sd = ip dd + xi, we may see 
straightforwardly that the phase dynamics for both ip dd 
and (p ss are closely related through the soft modes xi an d 
Xi in the two adjacent S layers. We simplify the equation 
of motion for ip ss by writing that ip' 



dd 



V>° 



Xl +X2- 



This substitution yields a familiar form of the equation 
of motion for both (p ss and (p dd . Here, we focus on the 
phase dynamics of ip ss given by 



Aire 2 



d 2 ip s 



d 2 <p s 



dt 2 



2 



jsd 

— — \sm(ip s 
e 



•EC- 1 ) 



. d 2 xi 

M dt 2 



C 



X 2 ) +sin(tp' 
9 2 X^ 



2J S 



Xi)] 



■ sw.(p 



(22) 



J; 7 



9x 2 



■ sin Xi 



= 



As indicated by Eq. ([22]) . the phase dynamics of v? ss 
depend strongly on the relative phases (i.e., xi an d X2) 
of the two condensates. As the phase Xi depends only on 
the intcrband Josephson effect, it is not influenced by the 
dynamics of the phase difference ip ss across the junction. 
Hence we may separate the phase dynamics decribed by 
Eq. (|22l) into two separate equations. Assuming that 



Xi 



= X2 : 

dx 2 



X, we write these equations as 
d 2 tp ss ( J sd 



cos x I sin if 



rq ^ — sm x = 



, (23) 
(24) 



d 2 x _ d^x 

dx 2 di 2 

where (£, x) and (£, x) are the dimensionless coordinates 
for the phases ip ss , and x^ respectively. Here, it is noted 
that, we may assume that ip ss = ip dd when xi — X2- This 
suggests that the two tunneling channels (i.e., tp ss = (p dd ) 
become equivalent. 



III. INTERBAND JOSEPHSON EFFECT 

In this section, we derive the equation of motion for the 
interband phase difference x of Eq. (|24p from the two- 
band Hamiltonian Htb.i of Eq. (TT]). It is noted that this 



Hamiltonian accounts for excitation of an i-soliton due 
to the interband Josephson effect. By following Leggettj 3 - 
we write the two-band Hamiltonian in the k-space rep- 
resentation as 



H 



TB,, 



i—s^d 



ni- 



OB,, 



E n 



yjjpair 

inter. £ 



■ E n 



rip 



(25) 



where the Hamiltonian T-i l OB g for the electronic band i is 



given by 



U 



OB,£ 



E 



Eh-C 



k c ka,£ C ktr,£ 



Ji J\ 



■Vu 2_^ c kt,^ c -k4,/ c -k'^,£ c ik;'t^ 



k,k' 



(26) 
Here, we introduce E so that the ground-state energy 
of ^\ Hq B 1 — E in the normal state vanishes for an 
arbitrary value of N s and N dl which denote the number 
of electrons in the s and <i-band, respectively. The inter- 
band pairing interaction between the electrons in s and d 
band is decribed by the Hamiltonian Titter 1 which may 
be written as 



inters 



-Vsd 



it if 
C kf,^ C -k^, 



C 1 C 1 

,c -k'V c k't,« 



(27) 



k,k',^j 



The term E ng denotes the ground-state energy of the 
system in the normal state. This energy E ng is fixed 
by the total number of electrons N = N s + N d in the 
system. The Hamiltonian T-Lf accounts for the effects 
of charge fluctuations from the equilibrium state. These 
charge fluctuations may arise from either the boundary 
conditions, various residual scattering processes, or the 
chemical potential fluctuations and lead to the charge 
imbalance between the s and d-band. The Hamiltonian 
"Kg 1 may be approximated 3 as 

Hf = lo [(N s - 7V S °) - (N d - N° d )\l = l(SN e ) 2 (28) 

when the deviation from the equilibrium is small. Here, 
7o = (pj 1 +Pd 1 )/8, Pi = m°kF/Tt 2 denotes the density of 
states for the i-band electrons at the Fermi surface, and 
N° denotes the number operator for the i-band electrons 
in the equilibrium state. 

To assess the dynamics of relative phase xe due to 
the interband Josephson effect, we rewrite the two-band 
Hamiltonain %tb,£ in terms of the familiar Ginzberg- 
Landau (GL) free energy. (See the appendix for a de- 
tailed discussion on the derivation of GL free energy for 
the two-gap superconuctors.) Here, we note that the tem- 
perature range is restricted to (T c — T)/T c <C 1, where 
T c is the superconducting transition temperature. Also, 
we considered the zero field limit here since the inter- 
band phase fluctuations do not carry magnetic flux. In 
the absence of magnetic field, we may write the one-band 
Hamiltonian H OB e in the form of the GL free energy 3 - 3 . 

Gob 1 as a function of pseudo-order parameter <&\ (i.e., 
Yl,iH- OB £ - E -> J2iGoB,i)- Here > % is tn e auxilary 



field representing an electron pair. For one spatial di- 
mensional case, we may write Q l OB « as 



'OB,. 



OB,e 

h GL 
,f£|fcJ|2 + 0L_| $ <|4 + 



1 

2m* 


dx 



(29) 



where af L and bf L are the coefhcients of the GL free 
energy and m* is the effective mass of the electron in 
the i-band. (See appendix.) When the order parame- 
ter §\ is written in terms of the modulus-phase variables 
[i.e., $^ = |$J|exp(i0|)], the interband pairing interac- 
tion Titter i yields the contribution 



frpair 



2.9. 



sd 



inter, £ 



9s9d 



\m\S>i\co S (e s t-6i), (30) 



where g t = (V ss Vdd - V? d )/Vu and g s d = V s d{V ss Vdd - 
Vsd) / (yssVdd) ■ Here, we assumed that the intraband in- 
teraction is much larger than the interband interaction 
(i.e., Vij < 7ji), With this assumption, we may approx- 
imate that |$|| ~ \/N a , which specifies that the intra- 
band interaction is independent of the coordinates (i.e., 
d\$>g\/dx sa 0). The charge- imbalance contribution %f 
in Eq. (J25I) may be expressed in a familiar form by using 
the number-phase uncertainty relationship of 



[sm, ( 



-4id 



w 



(31) 



where 6w is the kronecker delta. Noting that the Heisen- 
berg equation of motion for the phase difference xi — 
9f — Of is given by 



—jT =HXt>rlTB,e 



(32) 



we may write the charge imbalance contribution as 

or = — ( -^ \ ■ (33) 
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Combining the result of Eqs. (|2"5]l. (|50| and 
write the Gibbs free energy density Gtb z 



we may 
L~,ix3oB ~ 



Qng) + GinZr + Ci f° r tne two-gap superconductor in the 
superconducting state as 



Gtb,£ = / J 



a ? L m 



h GL 

+ ^-|^| 4 



I* 



' I 2 



2m* 



<w 



dx 



1 



dxt 
dt 



2g s 



\$l\\$j\ cos X t ■ (34) 
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Here, we assumed that the effects of external fields are 
weak inside each S layer and, thereby, their contributions 
in Gtb, i are neglected. 

We now use the Gibbs free energy of Eq. (fM)) and 
obtain the equation of motion for X£- By noting that 
there is no supercurrent Ji flowing anywhere inside the 
bulk superconductor, we write that 



Ji = 






m* dx 



= 0. 



(35) 



From the condition that Ji 
obtain that 



0, it is straightforward to 



m 



m 



+ 5 TT. 



(36) 



Here, 5 = for Jinter < (S-\ symmetry) and 8 = 1 

for Ji n ter > (S++ symmetry). We minimize the Gibbs 
free energy density with respect to the phase variation 
(i.e., dGTB,t/dxi = 0) and obtain the equation of motion 
for xi as 



d 2 x.e d 2 xe 



dx 2 



dP 



sin xi = , 



(37) 



where the dimensionless coordinates (x, t) are x = 



[45 srf «l* 



<'i - 



i d \<P s \ 2 )/g s g d \<£ s \\<P d \} 1/2 x and t 



t/(128j g sd \® s \\$ d \/g s g d ) 1/2 . It is straightforward to sec 
that Eq. (|3T|) may yield a number of soliton solutions, 
but, for simplicity, we consider a single-soliton solution 
of 



Xt(x,i) = 4 tan 



exp ± 



v t 



Vi 



(38) 



which is a kink-solution, known as the i-soliton. Here v 
is the speed of i-soliton. The kink-solution of Eq. (|38l) is 
identical to the functional form of the unperturbed fluxon 
solution. 24 However, the property of i-soliton differs from 
that of the fluxon. Unlike fluxon, an i-soliton does not 
carry a quantum of magnetic flux and cannot be driven 
by the Lorentz force due to the bias current J B as shown 
in Fig. [I] 

The equation of motion for xt and its solution of Eq. 
([551) indicates that the perturbation effects in the two- 
gap superconductor may lead the relative phase Xi to 
fluctuate from the phase-locked state of xi- These phase 
fluctuations, in turn, yield collective excitations, but they 
modify neither the ground state nor the one-particle ex- 
citation spectrum. However, if large phase fluctuations 
representing a 27r-phase texture can be stablized, then 
the one-particle excitation spectrum may become modi- 
fied. 



IV. EFFECTS OF i-SOLITION ON PHASE 
DYNAMICS 

We now examine the effects of large fluctuations in the 
relative phase of the s and d condensates on the fluxon 
dynamics. When the amplitude of these phase fluctua- 
tions grows to the nonlinear region and becomes stablized 
as suggested by Tanaka^ excitation of an i-soliton can 
change the amplitude of critical current density. The 
modulation of critical current between two adjacent S 
layers may be considered as an effective potential for a 
fluxon. The effects of an i-soliton excitation on the phase 
dynamics of the LJJ may be described by Eq. (|2~3")) . We 
may simplify Eq. &j§ by substituting Eq. ([3"g]). The 




FIG. 2. The amplitude modulation of the Josephson current 
density J c (x)/J° due to excitation of static i-soliton repre- 
senting 27r-phase texture with its center located at x = 
is illustrated for the interband Josephson current density of 
jsd/jss _ 1Q ( solid line ) an d _o.io (dashed line). Here, J° 
is the critical current density of the homogenous LJJ in the 
absence of the interband Josephson effect. 



equation of motion for the fiuxon which accounts for the 
interband Josephson effect in each S layer is given as 



dx 2 



d 2 <p 

dt 2 



Jc . 



0, 



(39) 



where <p = ip ss , J c /J° = 1 + (J sd /J ss )[l - 2sech 2 (a x - 
/3 t)] and J° is the critical current in the absence of 
the interband Josephson effect (i.e., J sd = 0). Here 
«„ = [J^Xl/il-vDJ^dsd] 1 / 2 &nd(3 = v [J sd fi 2 d e/(l- 
v 2 )J ss d s d] 1 / 2 . Equation (|39)) indicates that a i-soliton, 
representing a moving 2tt- phase texture of Eq. (|38p , leads 
to both spatial and temporal dependent modulation of 
the critical current. 

The effects of spatial and temporal variation of J c /J° 
on the fiuxon depend on the shape of this modulation. 
In Fig. [21 we neglect the temporal modulation and plot 
the amplitude of the critical current (J c /J°) versus the 
dimensionless spatial coordinate x for J sd /J°=0.1 (solid 
line) and -0.1 (dashed line) to illustrate the effects of a 

single z-soliton excitation in the S++ and S^ symmetry 

superconductor, respectively. Here, we set a = 1 for 
definiteness. The curves show that the shape of critical 
current modulation depends on the symmetry of order 
parameter (i.e., S++ versus S-\ — ). However, as we discuss 
in Sec. V, when J sd / J ss -c 1, the symmetry of the order 
parameter does not affect the fiuxon motion significantly. 

The critical current modulation induced by the inter- 
band Josephson effect has two main effects. First, the 
shape of fiuxon may become deformed. However, for a 
small modulation considered in the present work, this ef- 
fect is negligible. Second, the speed of fiuxon becomes 
modified since the critical current modulation behaves 
as an effective potential. In the region of critical current 
modulation, the fiuxon speed may become significantly 



changed from a uniform value. These changes imply that 
the EM waves can be emitted by a moving fiuxon as it 
decelerates. 

We use the perturbation method to examine the ef- 
fects of critical current modulation on the emission of 
EM waves. To obtain physical insight, we carry out the 
calculation in the rest frame of the fiuxon (i.e., a reference 
frame which is moving with the speed of the unperturbed 
fiuxon) as described^ by Fogel and coworkers. To this 
end, we perform the Lorentz transformation of 



F = 



t — vx 



vT 



and 



x — vt 



VT 



(40) 



where v is the speed of the unperturbed fiuxon. With this 
transformation, we rewrite the sine-Gordon equation of 
Eq. ([29J as 



d 2 <p 
dx' 2 

1 - 



dt' 2 



Jc 



■'c 

{J sd /J ss )[l 



smip = , (41) 

2sech 2 (a' x' + f3' i')}, 



where J c / J° 

a' a = (a - (3 v)/VT^ 2 ~ and j3' = (a v - /? )/vT^A 
By considering the case of weak interband Josephson 
effect (i.e. J sd /J ss < 1), we assume that a solution to 
Eq. (14"TT) may be written as 



jsd 
ip(x,t) W (p (x) + —ipi(x,t) 



(42) 



Here, for convenience, we make the following changes in 
the notation: (x',t') — ¥ (x,t). A solution to the unper- 
turbed sine-Gordon equation of 



d 2 ip d 2 ip 



dx 2 



dt 2 



sin (p = 



(43) 



is given by ip (x) = 4tan~ 1 [exp(x)]. Following Fogel et 
al.,— we may write the correction term ipi due to the 
critical current modulation in the most appropriate basis 
by separating the spatial and temporal dependence as 



(fi(x,t) = f(x)e~ 



(44) 



The separation of variables for the perturbative contri- 
bution cpi in the rest frame of the fiuxon (i.e., v = 0) 
leads to the eigenvalue equation for f(x) as 



— + (1 
dx 2 



2sech x) 



f{x) = u 2 f(x) 



(45) 



The eigenvalue problem of Eq. (|45|) yields one bound 
state with to = u>b = and a continuum of scattering 



states with to = us 
malized eigenfunctions are 



1 + k . The corresponding nor- 



m 



' fb{Xy - 

J{k,x) = 



2sechx , 

"-t-itanhic inx 



with u = W& 
with <jj = lu k 



(46) 



Here we use the subscripts b and k to denote the bound 
state and continuum of scattering state k, respectively. 
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The bound state fb(x) is associated with the Goldstone 
translation mode of the fluxon, while the continuum 
eigenfunctions /(k, a;) represent the radiation modes. 
The eigenfunction of Eq. (|4"6")l indicates that the first 
order correction (pi(x,t) due to the critical current mod- 
ulation may be separated into two parts as 



<Pl(x,t) = tptrans(x,t) + (fi rad (x,t) 



(47) 



Here, tptrans and if ra d represent the bound state and con- 
tinuum eigenstate contribution, respectively. 

The bound state contribution iftrans(x, t) may be writ- 
ten as 



iptrans{x,t) = -(j) b (t) f b (x) 



(48) 



The amplitude (f>b(t) of the bound state is determined 
straightforwardly from the equation of 



d 2 db(t) . f°° , , „ , o* -, sinha; 
-4/ dx(l - 2sech 2 £ ) — - T - 



dt 2 



sech x 



(49) 



where £ = a' x + f3' t. The solution to Eq. (l4T)l) may be 
obtained as 



hit) = -^S I i 



ft 2 



dx- 



sech x 



2P' t e 2a' x I I 



(50) 



We note that 4>b{t) may be used to evaluate the trans- 
lation mode contribution Lp trans {x 1 t). This contribution 
has no effects on the motion of the fluxon center. 

The continuum eigenstate contribution, representing 
the radiation modes, is given by 



<Prad(x,t) 



dn (j>(K,t)f(K, x). 



(51) 



The amplitude <P(k, t) is determined from 
d 2 cf)(K, t) 



dt 2 



+ (l + K 2 )<j)(K,t) = Q(n,t) , (52) 



where 



/°° sinliT 
dxf*( K ,x){l-2sech 2 £, ) =-. (53) 
-co cosh X 

In obtaining Eq. (|52j) . we used the orthonormality con- 
dition for the eigenfunctions (i.e., J dxf*{n' ,x)f(n, x) — 
5(k — k')). The contribution to the radiation mode of tpi 
may be estimated by solving Eq. (|52[) . For a single mod- 
ulation of the critical current density as shown in Fig. 
[H we may obtain a solution to Eq. (|5"2")l more easily by 
using the relation 



sech £ 



dk irk 



2ir sinh 



Trk 



sKo 



(54) 



which is the Fourier representation of the critical current 
variation. Using this substitution, we rewrite Q(ft, t) by 



integrating the right hand side of Eq. (f53|) over x and 
obtain 



Q(M) 



, (1 + K 2 )sech- 

r.rch 7rfe7 " , 1 

dfc 2 , [1 + k 2 - (ka) 2 ]ke ik ^ ,(55) 

2sinhf J A 

where r) K = a' Q — (n/k). The solution 0(k,£) may be 
written as 



(M) = 



dw Q(k, w) 

2^ (1 + K 2 ) - OJ 2 



(56) 



where Q(k, w) = J dt'Q(K,t')exp(—iuit'). It is straight- 
forward to evaluate the integration over t' and uj, and 
write the solution <j)(n, t) as 



(«,*) = 



^2^(1 + ^2) 
sech^ 



K7T 

sech — 
2 



difc- 



(fc^) 2 



sinh^ 1 + k 2 - (fc/^) : 



jfceifc^t 



(57) 



However, as indicated in Eq. (|5"Tj) . we need to integrate 
over the continuum variable k to compute <p ra d(x, t). 
This may be evaluated by using the contour integration 
method. The location of the poles for the contour integral 
is shown schematically in Fig. [3] The radiation contri- 
bution ip ra d(x,t) of Eq. (f5i"|) indicates that all poles are 
simple, and the residue of each pole may be evaluated 
separately. The loca tion of the poles are the following: 
z = +i, z x = +i\/l - (k/3' ) 2 , z° n = +i(2n + 1), and 
2„ = ±ka' +i(2n+l), where n = 0, 1, 2, •••. The residues 
of the pole structure shown in Fig. [3] yield two types 
of contribution: i) exponentially localized contribution 
around the fluxon center and ii) linear traveling wave con- 
tribution. Hence, we may decompose the radiation mode 
ip r ad into the exponentially localized (frad) ana - traveling 
wave [VraD contributions: <p rad = <p%» + <^ e . We 
note that the exponentially localized contribution ^Vad 
does not produce a true radiative correction. Only the 
traveling wave </?™ a J e gives rise to a true radiative contri- 
bution. 

We now focus on the poles that give rise to the trav- 
eling wave contribution to the radiation correction. As 
indicated by the Fourier components of the critical cur- 
rent modulation described by the (1 — 2sech £ ) factor in 
Eq. (fS"3")) . the condition for the traveling wave radiative 
contribution depends on k. For fixed k < (l//3' ), the pole 
at z\ lies on the imaginary axis as shown in Fig. [3J How- 
ever, as the fluxon speed v increases, the pole z\ moves 
down the imaginary axis. At the critical value v = v t h, 
the pole z\ lies at the complex plane and it becomes real 
when v > Vth ■ The changes in the radiation contribution 
in Eq. ([51]) from this pole may be easily identified by the 
contour integration since it is not exponentially localized 
around the fluxon center but oscillates with x. The oscil- 
latory contribution only arises when v > v t h- This leads 
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FIG. 3. The pole structure for the radiation contribution 
of Eq. (|51|) to (pi (k, t) for a fixed k is shown schematically. 
The solid circles represent the poles yielding the exponentially 
localized contribution. The open circles represent the changes 
in the location of Z\ as the fluxon velocity v changes from 
V < vth (on the imaginary axis) to v > Vth (on the real axis). 
The shift direction for the pole z\ is indicated by the arrows. 



to the radiative contribution of 



^ctom) = - 



P'o 2 



D dk Ttk 
2r] sinh^ v ' ' 

a ik(p' t-\-r)x) p —ik(/3 t—rfx) 



cosh 



Trkr/- 



cosh 



■7rkr}-\. 



(58) 



where r\ = \J(k(3' ) 2 — 1, and r]± = a' =p (r]/k). Equation 
(jS"8f indicates that, for a fixed k, this radiation correction 
is the superposition of two linear traveling waves with 
different amplitudes. The two waves travel in opposite 
directions. The threshold velocity Vth for the fluxon is 
given by 



Vth 



[4^ 
\ 2 



1/2 



(59) 



The dependence of Vth ° n the z-soliton velocity v indi- 
cates that, for the case of static spatial variation of the 
phase (i.e., v = 0), EM radiation may be emitted by the 
fluxon whenever it passes through a region where the crit- 
ical current is affected by the interband Josephson effect. 
Hence, when an array of static i-solitons are excited to 
yield a spatially periodic modulation of the critical cur- 
rent density, the threshold velocity v t h becomes finite. 
This radiative threshold is similar to that found in ear- 
lier studies 2£r— 



V. CURRENT- VOLTAGE CHARACTERISTICS 

Emission of EM radiation by a moving fluxon, indi- 
cated by the linear traveling wave contribution to the 
radiative correction in tpi(x, t) as discussed in Sec. IV, 



reflects the changes in the fluxon dynamics. We now ex- 
amine the effects of a single i-soliton excitation on fluxon 
dynamics by considering a static phase texture described 
by Eq. ([M]) with v — (i.e., j3 = 0). Here, we include 
the perturbative effects of bias current and dissipation to 
examine a realistic tunnel junction. The bias current J B 
acts as a driving force for the fluxon, while the two dis- 
sipation terms, Ti(dip/dt) and T 2 {d 3 tp / did 2 x) , account 
for the interaction between the fluxon and dissipative en- 
vironment. Using the perturbed sine-Gordon equation 
with the critical current density modulation, we deter- 
mine the fluxon trajectories and estimate the effects of 
the bias current and dissipation (Ti and 1^). These per- 
turbation contributions as well as the interband Joseph- 
son effect can modify the current-voltage curve. We start 
with the perturbed sine-Gordon equation of 



d\p_ 

dx 2 



dt 2 



sirup = ^(ifjXjf) 



(60) 



where the perturbation term T((p, x, i) is given by 



jsd 

T{<p,x,fj = — [1 



jB 

sech (a x)] simp — 

J r 



dip 



r 5 



d z tp 
dtd 2 x 



(61) 



Here, we assume that each perturbation term in J- is 
small and does not change the shape of the fluxon in 
the leading order. The main effect of the first term of 
T in Eq. (|6ip is to provide a potential for a moving 
unperturbed fluxon of 



<p (x, i) = 4tan 1 (e 1 ') , 



where 



C(x,F) = ± 



x - fo v(t)dtf - Xoffl 



v/I - ~& 



(*) 



(62) 



(63) 



Here, the fluxon speed v(t) accounts for the time depen- 
dence induced by the critical current modulation. We 
now examine the trajectories of the fluxon in an LJJ, 
and the effects of an i-soliton excitation on the I-V curve 
in the low-voltage regime. 

By following McLaughlin and Scott ^ we examine the 
fluxon trajectories by computing the fluxon speed and the 
corresponding position. For this purpose, we assume that 
the fluxon approaches the region of a large stable varia- 
tion of phase difference (i.e., 2-7r-phase texture) between 
the s and d condensates from x = — oo with increasing t. 
We write Eq. (JHOJ) as two first order differential equations 
describing the velocity (v) and the position (x ) of the 
fluxon, respectively, as 

dv 1 - v 2 f°° 

=F — -. — / dx T((p ,x,t) sech£ , (64) 



dt 



dt 4 



dx F(ip ,x,t) £sech£, (65) 



where £ = ((x,t) is defined in Eq. (J63I) . The pertur- 
bation term T modifies the speed v of the unperturbed 
wave form of Eq. (|62j) to depend on time t. By perform- 
ing the integration, we rewrite Eqs. (|64[) and (|65[) . which 
respectively account for the fluxon speed w and position 
X, as 



-= = ± (1 - V 2 )? - Tivfl - v 2 ) - -T 2 v 

dt 4 JS v ; v 7 3 



+ (1- W 2 )' 



,/ 



sd r 00 



dyT. 



sinhy 
cosh 3 y 



(66) 



and 

dX i 

-£ =v+ ~ 

where Y± = 
sech 2 a (vT" 



v 6 J s 



J s 



1-2 



*T_^]. ,07, 
cosh y 



T±(y,t,X) = sech a D (Vl - v 2 y + X) ± 
^2y - X) and X = f*v{t')dt' + x (t). 
These two equations describe the fluxon trajectories in 
the (v, X) phase plane. We numerically integrate Eqs. 
(IBS]) and (J6"T|) to estimate the fluxon trajectories. The 
uniform fluxon speed v = v^o far away from the region 
of critical current modulation is given by 



dv _ _ i _ 7rjB ( 



v 2 )* -riu(i-u 2 ) 



i 

—Tov. 
3 



(68) 



The power-balance velocity Voo may be estimated by set- 
ting dv/di = 0. The fluxon speed Voq is obtained by 
solving the following cubic equation: 



i=0 



= 



(69) 



where z = (1 - v^) 1 / 2 , a 3 = [^ 2 (J S ) 2 /16(J C °) 2 ] + T 2 , 
a 2 = (2rir 2 /3) - T 2 , en = (r 2 /9) - (2rir 2 /3) and a Q = 
— T 2 ,/!}. We note that the solution is bounded by the 
condition that < Uoo < 1 since Uoo is given in units of 
Swihart velocity. 

We now estimate the effects of bias current J B / J° on 
the fluxon trajectories in the (v, X) phase plane. In Fig. 
|4j the fluxon trajectories obtained by numerically inte- 
grating Eqs. fl66|) and (J6TJ) are plotted for J B /J° = 0.01, 
0.02, and 0.04 (from left to right) to illustrate the posi- 
tion X of the fluxon as a function of velocity v. Here, 
we set the dissipation parameters Ti = T 2 =0.1. The 
solid (dashed) curves represent the interband Joseph- 
son current density J sd /J ss = 0.1 (-0.1). Here, Uqo is 
the uniform initial speed of the fluxon at a position far 
away from the region of 27r-phase texture which is cen- 
tered at X = 0. The value of -u^ depends on J B / J°. 
The curve for each J B / J° shows that, as the fluxon ap- 
proaches X = 0, the fluxon speed deviates from a straight 
vertical dotted line representing a uniform speed in the 
absence of the critical current modulation. The curves 
also show that when the bias current is small (i.e., see, 
for example, J B / J° — 0.010) the fluxon becomes almost 
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FIG. 4. Fluxon trajectories in the (v, X) phase plane for 
Ti = T2 =0.1 are plotted to illustrate their dependence on 
the strength of both the bias current J B /J° and interband 
Josephson curent J sd /J SB . Three solid (dashed) curves, from 
left to right, correspond to J B /J° = 0.010, 0.020 and 0.040 
for J sd /J ss = 0.10 (-0.10). The vertical dotted lines represent 
the uniform fluxon speed in the absence of critical current 
modulation. 



pinned since the pinning effects of the critical current 
modulation reduce the fluxon speed close to zero. For 
J B /J° < 0.100, the fluxon is pinned: t) = 0atl« 6.5 
for J sd /J ss = 0.10 and X w -6.5 for J sd /J ss = -0.10. 
However, for J B / ' J° > 0.010, the fluxon is not pinned 
since the driving force due to the bias current is larger 
than the pinning force. The fluxon, approaching from 
X = — cx3, undergoes a notable change in its speed in the 
region of the critical current modulation. The difference 
between the maximum and minimum value of the fluxon 
speed becomes reduced with the increasing bias current 
density. This is due to the fact that the importance of 
the pinning effect decreases with increasing value of v^, 
(i.e., increasing J B / J°). 

As its speed decreases, the fluxon radiates EM waves 
in the region of critical current modulation. A variation 
in the fluxon sped may be seen easily in Fig. [J] as a devi- 
ation of the fluxon trajectories from the vertical dotted 
lines. For the solid lines representing the S++ symmetry 
(i.e., J sd /J ss > 0), the fluxon speed v > Voa for X < 0, 
but v < Woo for X > 0. On the other hand, for the dashed 

lines representing the S^ symmetry (i.e., J sd /J ss < 0), 

the fluxon speed v < v^ for X < 0, but v > Woo for 
X > 0. The deviation from the uniform fluxon speed 
Voo increase with the interband Josephson current den- 
sity \J sd /J ss \. This non-uniform fluxon speed due to 
the critical current modulation suggests that the fluxon 
absorbs energy from the interband Josephson current to 
increase the speed from v 00l but it returns the energy 
back to the junction by radiating EM waves, resulting 
the decrease in the speed. 

As suggested by the fluxon trajectories in the (v, X) 
phase plane, the bias current density must be larger than 
the threshold value Jf h in order for the fluxon to pass 
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FIG. 5. The current- voltage curves for J sd / J sa = 0.10 (dotted 
line) and 0.20 (solid line) illustrate the dependence of the 
threshold bias current (i.e., the value of J B / J° for V — 0) 
on the interband Josephson current. Here, the dissipation 
parameters are Fi = 0.10 and T2 = 0.10. 




FIG. 6. The threshold bias current density jf K i ' J° is plot- 
ted as a function of J sd /J ss for Ti = 0.05 (solid line), 0.10 
(dashed line), and 0.15 (dot-double dashed line) to illustrate 
the minimum bias current needed to overcome the pinning 
effect of the critical current modulation. Here, F2 = 0.10. 



through the region of critical current modulation. In Fig. 
U we plot the I-V curve for J sd /J ss = 0.10 (dotted line) 
and 0.20 (solid line) to illustrate the dependence of the 
threshold bias current jf h on the interband Josephson 
current (J sd /J ss ). The fluxon speed, according to the 
AC Josephson effect (l/i)(dip/dt) — V, is proportional 
to the voltage V across the junction. Here, we set the dis- 
sipation parameters as Ti = T2 = 0.1 for concreteness. 
The value of the threshold bias current Jf h for a fixed 
jsdjjss ma y kg eg ti m ated easily since the voltage across 
the LJJ does not appear until the bias current reaches the 
threshold value. The curves show that the threshold bias 
current increases with J sd / J ss . As the critical current 
modulation plays the role of an effective potential for the 
fluxon, a larger J B / J° is needed to overcome the pinning 
effect as J sd /J ss increases. Hence the threshold bias cur- 
rent is similar to the minimum current density needed to 
overcome the pinning force. In Fig. [SJ for J B > jf h , 
the I-V curves show that the voltage increases steadily 
with increasing J B because the boundary effects are not 
included in this work. The zero-field step resonances due 
the boundary current^ (i.e., a rapidly vanishing critical 
current near the edges of the junction) are expected to 
be present along with a smooth increase shown in Fig. [S] 
when the boundary effects are included. 

In Fig. H we plot J B / ' J° as a function of J sd /J ss 
for Ti = 0.05 (solid line), 0.10 (dashed line) and 0.15 
(dot-double dashed line) to illustrate the dependence of 
the threshold bias current on the interband Josephson 
effect. We set T2 = 0.10 for concreteness, but the curves 
are independent of the strength of the parameter T2- The 
curves show that Jf h increases linearly with the interband 
Josephson current J sd /J ss . As indicated by the fluxon 
trajectories in the (v, X) phase plane, the J B h j ' J° versus 

jsd J JSS ^^gg for both jsd I JSS > and jsd I JSS < Q 

are identical. This indicates that the symmetry of or- 



der parameter (i.e., S ++ versus S-\ ) may not be dis- 
tinguished directly from the TV curves. The increase 
in the threshold bias current with increasing interband 
Josephson current implies that, as the critical current 
modulation increases, a greater strength of driving force 
must be provided by the bias current to allow the fluxon 
to pass through the region of a static 27r-phase texture. 



VI. SUMMARY AND CONCLUSION 

In summary, we have investigated the effects of inter- 
band Josephson current on the fluxon dynamics of an 
LJJ with two-gap superconductors. Due to the multi- 
component nature of the superconducting order param- 
eter, a 2-7r-phase texture may be present in each S layer. 
Accounting for the charge imbalance between the two 
electronic bands, say s and d bands, we found that the 
dynamics of the relative phase of the s and d condensates 
may be described by the equation of motion which is sim- 
ilar to that for the fluxon. However, the kink-solution 
for the relative phase (i.e., i-soliton) differs from that for 
the fluxon. Unlike fluxon, an i-soliton does not carry 
a magnetic flux quantum but leads to both spatial and 
temporal modulation of the critical current density. 

The critical current density modulation, induced by 
excitation of an z-soliton representing a 27r-phase tex- 
ture, behaves as an effective pinning potential for fluxon 
and modifies the fluxon trajectories. These trajectories 
in the velocity-position phase plane are similar to those 
for fluxon moving in the presence of a single microre- 
sistor. Moreover, similar to an array of microresistors 
in the insulator layer of a LJJ, a periodic modulation 
of the critical current may be created when a spatially 
periodic array of i-solitons is excited. This critical cur- 
rent modulation can change the fluxon speed when the 
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bias current is applied across the junction in a dissipative 
environment. In the region far away from the critical cur- 
rent modulation, the fluxon motion is uniform. However, 
in the region of ciritical current modulation, the fluxon 
speed can either decrease or increase from the uniform 
speed as the fluxon either receives or returns energy to 
the junction. The decrease in the fluxon speed due to 
excitation of an i-soliton results in emission of EM radi- 
ation (i.e., quasi-linear plasma wave). 

The changes in the fluxon speed due to critical cur- 
rent modulation may be reflected in the I-V character- 
istics. When the bias current is less than the critical 
value, the fluxon becomes trapped and its speed will be 
reduced to zero, resulting zero voltage across the junc- 
tion. This is similar to a fluxon trapped by a microre- 
sistor in the I layer. The result suggests an interesting 
possibility that a double- well potential for a fluxon, which 
is a necessary condition for a Josephson vortex quantum 
bit, may be created by exciting two i-solitons in each 
S layer of the LJJ, rather than implanting two closely 
spaced microresistors^ Also, the I-V curves reveal the 
dependence of fluxon motion on the bias current. The 
threshold bias current needed for the fluxon to overcome 
the critical current modulation appears as a discontinu- 
ous jump in the I-V curve at V=0. Hence, this discon- 
tinuity may serve as a way to verify the excitation of an 
i-soliton. Since the size of this discontinuity depends on 
the interband Josephson current, the I-V curves may be 
used to estimate the strength of the interband Josephson 
effect. However, the I-V curves for both the S++ and 
5-| symmetry superconductors are expected to be sim- 
ilar since the threshold bias current does not depend on 
the symmetry of the order parameter. 

The present work indicates that the changes in the 
fluxon dynamics due to excitation of an i-soliton in a 
single LJJ may become amplified in a multiple LJJ stack 
where the collective motion of fluxon is known to arise. 36 
Furthermore, as recent studies of three-band supercon- 
ductors indicate, excitation of j-solitons 3 -^ which breab 3 ^ 
the time reversal symmetry is possible when the inter- 
band Josephson couplings are frustrated^ The fluxon 
dynamics of a LJJ with three band superconductors may 
lead to a more interesting junction property than that 
described in the present work. 
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APPENDIX: 

DERIVATION OF GINSBERG-LANDAU ENERGY 
FOR A TWO-GAP SUPERCONDUCTOR 



For completeness, we derive the Ginsberg-Landau free 
energy from the two-gap Hamiltonian T-Ltb of Eq. ([1]) by 
writing the partition function Z as 



V[c s ,c s ]V{c a ,c d ]e 



where the action S is given by 



d „di-S[c*,c°,c d ,c d 



[P r 

JO L ; : ,, _ 



(70) 



(71) 



d T = d/dr and c l ka (c kcr ) is the Grassmann variable rep- 
resenting the creation (annihilation) operator for the i- 
band electron in the (k<r) state. Here, we suppress the S 
layer index £, for simplicity. We proceed by introducing 
Nambu notation of C k = {Sj s *-c l _j s i ) and by writing the 
pair fields A\ and A\. as 



-4k = Cl T +Cl = Cktc-ki 



(72) 
(73) 



where t± = (n ± iT-^jl and Tj are Pauli matrices. The 
action S of Eq. ([Til may be rewritten as 



Ja i,k L j,k> 



(74) 



We introduce the Hubbard-Stratonovich transforma- 
tion to map the interacting system to non-interacting 
fcrmions moving in an Hubbard-Stratonovich field $ (i.e., 
auxilary field) representing electron pairing. We rewrite 
the partition function of Eq. (fTDj) as 



Z= 2>[C,C]2?[*,*]e 

where the action S is given by 

S = I dr\Y,Ci{d r I + A^)Ci 

Jo L .- 1, 



-s[c,c,*,*] 



(75) 



W^c^k'-A-^lk' 

k.k' 



(76) 



Here, we note that $ k = (<F£, $£), _4 k = (A£,A£), and 
V denotes the pairing interaction matrix 



V = 



v ss v sd 

V s d Vdd 



(77) 



We shift the $ field (i.e., <& k ->• $ k + A^V and $ k 
^k + V^k) and obtain 

r0 



S= f dr £C k (3r/ + e k T 3 )C k 
Jo L i,k 

+ Yl ( *kTT$k' + A$k' + $k-4 k ' 



k : k' 



(78) 



13 



We integrate out the fermion variables C and C by using 
the Grassmann integrals and obtain 



S 



& 



k.k' 



y-^kik 



9sd 

g s gd 



(K*i>+nK>) 



-TrlnG*; 1 -TrlnG^ 1 



(79) 



where g t = (V ss V dd - V^ d )/V u and g sd = V sd {V ss V dd - 
V* d )/{V ss V dd ). We note that the first term in Eq. ^ 
may be considered as the pesudo-order parameter contri- 
bution Sgap to the action. The interacting Green function 
d for the i-band electron is given by 



g: 



(G°y 



(80) 



where the non-interacting Green function G° is given by 



{G°r i = 



Or 







d r - el 

and the pair interaction £j is given by 

$ 



*t 







(81) 



(82) 



We may simplify Eq. 
Tr In G7 1 as 



inGz^HGtr 1 -^ 



79]) by expanding the terms 



<G°x z y 



(83) 



In the expansion, all odd order terms vanish when the 
trace of InG" 1 is taken. Assuming that the pair field $^ 
is uniform, we may write the fourth-order contribution 
in the expansion of In G~ as 



d 3 k 



7C(3)a 



2|$ 4 



K) 2 ? 



4n 2 T 2 



I* 



(84) 



by making the high density approximation. Here, cj„ = 
(2v+l)ir/f3 with v = 0, ±1, ±2, ••• denotes the Matsubara 
frequency, and £(x) is the zeta- function The second order 
contribution to the expansion yields 



Tr(G°S,) 2 



E 



l*U<l 2 



/? ,lTk' ( ^~ e k)(^ + e k') 



(85) 



We may simplify the calculation by introducing new mo- 
mentum variables: p = (k + k')/2 and q = k k'. We 



make a Taylor expansion in q up to second order, in the 
high-density limit, and obtain 



Tr(G°E, 



^Ei^i 2 E 



1 



(^) 2 



+ 



2 9 

pZqZ 



12(m°) 



oY2 



4fe) 2 



.(^ + (4) 2 ) 2 ' (^ + (^) 2 ) 3 JI' (86) 

Here, we used the dispersion relation of e^ = k 2 /2m°, 
where m° is the band mass of the electron. Also we 
used the relation piPj — > Sijp 2 /3 when the angular av- 
erage is performed. In Eq. (|86|) . the first term yields 
the quadratic term (i.e., |<I> 1 | 2 ) in the GL free energy, 
while the second term of leads to the gradient part of the 
quadratic term (i.e., |V<I> 2 | 2 ). Combining these terms, 
we obtain 



^Tr(G°S,) 2 



^In^Vf 



7((3)k 2 F (3 2 Pl 
24(rn°) 2 7r 2 



d 3 r|V$ ,; (r)| 2 



(87) 



where A is the cut-off energy for the boson which medi- 
ates pairing interaction and 7 sw 0.5772. Here, we ap- 
proximated that p 2 fa kp. We now combine the result 
and write the coeffcients to the GL free energy expansion 

for gp B = Q l OB - 

ducting state as 



Q % for one dimension in the supercon- 



1 



'OB 



h GL 

«f L |$T + ^-|^| 4 -r — 



d$> 1 



d.r 



where the coefficients of the expansion are 



,GL 



W L = 



1 , 2/3Ae" 7 
pihi , 

9i 7T 

Kim 2 p l 



16tt 2 
and the effective mass m* is 



(89) 
(90) 

(91) 



The GL free energy of Eq. (1881) may be used to esti- 
mate the interband phase dynamics by noting that the 
Hubbard-Strotonovich field may be expressed in terms 
of modulus-phase variables as $i = |$ 4 | exp(i#*), as dis- 
cussed in Sec. III. 



o\2—2 



24(m°) 2 7r 
7C(3)fc|p 4 ' 
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